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We investigate the lifetime of macroscopic entanglement under the influence of decoherence. For 
GHZ-type superposition states we find that the lifetime decreases with the size of the system (i.e. 
the number of independent degrees of freedom) and the effective number of subsystems that remain 
entangled decreases with time. For a class of other states (e.g. cluster states), however, we show 
that the lifetime of entanglement is independent of the size of the system. 
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The question whether or not entanglement - a genuine 
feature of quantum mechanics - can persist in a macro- 
scopic (i.e. "classical") world, has entertained quantum 
physicists since Erwin Schrodinger introduced his notori- 
ous gedanken experiment known as "Schrodinger's cat" 
PjQ. While entangled states of microscopic matter, such 
as a few atoms or ions in a trap, can nowadays be pre- 
pared in the laboratory it is often argued that for 
a large number N of particles this task would become 
exceedingly difficult, since the effective decoherence rate 
would grow linearly with the size of the system N . 

In its simplest version, the argument is based on the 
observed evolution of superposition states of the form 
\GHZ) = l/^/2{\{))®^ + \\)®^), also called Greenberger- 
Hornc-Zcilinger (GHZ) states, of a system of N spins 
or qubits interacting with uncontrollable degrees of free- 
dom of the environment, described e.g. by a heatbath. 
The rate at which this state decoheres scales indeed as 
kN , where k is the decoherence rate of a single qubit. 
While this observation is correct, it is not clear whether 
the scaling of the decoherence rate with the size of the 
system iV is a special property of GHZ states or a gen- 
eral feature of all multiparticle entangled states. This 
means that the conclusions that are usually drawn from 
this observation, namely that macroscopic entanglement, 
i.e. entanglement between a macroscopic number of par- 
ticles, necessarily becomes exponentially fragile with N , 
are questionable and will indeed be refuted in this paper. 

To this aim, we investigate the effect of decoherence 
on the entanglement properties of a class of multiparticle 
entangled states. We consider the lifetime of entangle- 
ment between a variable number of subsystems of the sys- 
tem. Specifically, we consider both the time after which 
the distillable entanglement of the state vanishes and the 
time after which the state becomes separable. These life- 
times are, in general, finite. We determine the scaling 
behavior of these lifetimes with N . For GHZ states we 
find that an increasing decoherence rate indeed results in 
a lifetime of distillable A^-party entanglement 0, that 
decreases with N . For many other states (e.g. cluster 
states Q), in contrast, we show that the lifetime of gen- 
uine multiparty entanglement is independent of the size 
of the system N . In particular, we find that the lifetime 
of any state that belongs to the class of graph states 
\4>)g |y| -which contains the GHZ and cluster state as 



particular cases- is bounded from below by a quantity 
which depends on the maximum degree of the associ- 
ated "interaction" graph G, but is independent of N . 
This implies that genuine multi-particle entanglement of 
a macroscopic number of particles is possible and can 
persist for timescales that are independent of the size of 
the system. 

When describing multi-particle entanglement, a cen- 
tral notion is the partitioning of a system of N particles 
into M < N groups. Each group may consist of several 
particles, which are then considered as a single subsystem 
with a higher-dimensional state space. If we associate a 
specific spatial distribution with the particles, as in the 
case of spins on a lattice, partitionings may be chosen 
that correspond to a re-scaling of the size of the subsys- 
tems, as it is used in statistical physics. In the present 
context, we will be interested in the behavior of the distil- 
lable entanglement under a coarsening of the partitions, 
i.e. under a re-scaling of the size of the subsystem that 
is under local control of a single party. We will con- 
sider the case where the entire lattice of N qubits is in 
a graph state and compare in particular the two cases 
of GHZ states and cluster states. On the one hand, this 
approach allows us to determine the effective size of the 
system, i.e. the number of subsystems which are still en- 
tangled after a certain time. On the other hand, we can 
investigate the behavior of entanglement under re-scaling 
in the asymptotic limit N ~* oo. We find that the life- 
time of the distillable entanglement of the cluster state is 
largely independent of the size of the partitions, and thus 
the same on all scales. For the GHZ state, in contrast, we 
find that the distillable entanglement vanishes after and 
arbitrary short time on all scales, as long as we consider 
partitions of finite size. However, if we allow the sizes of 
the partitions to become macroscopic themselves (in the 
sense that the A'^ qubits are divided into a fixed number 
of M cells whose size N/M grows to infinity as iV ^ oo) 
then the lifetime of this AZ-party distillable entanglement 
(between M cells of macroscopic size) becomes finite and 
scales to leading order as 1/(kM). 

Throughout this paper, we will mainly use a decoher- 
ence model corresponding to individual coupling of par- 
ticles to a thermal bath in the large T limit, described (in 
quantum information language) by individual depolariz- 
ing channels. Although we are aware of the restricted 
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applicability of such a decoherence model, we emphasize 
that any model with independent couplings of particles 
to the environment (which appears naturally in the case 
where the N particles are spatially separated and thus 
interact with independent environments) can be mapped 
(by introducing additional noise) to such a depolarizing 
channel. This model — in contrast to the widely used de- 
phasing channel — is basis independent, i.e. invariant un- 
der local unitary operations. In addition, in many cases 
it turns out that the scaling behavior of the lifetime of 
entanglement with the size of the system is independent 
of the specific decoherence model. 

We consider depolarizing channels with noise param- 
eter p = p{t) = e^'^*, where k is a decay constant de- 
termined by the strength of the coupling to the environ- 
ment and t is the interaction time. The channel acting 
on particle k is described by the completely positive map 
(CPM) 



£kP = p{t)p- 



l-p{t) 



{k) (k) 



(1) 



where aj are the Pauli matrices with cto = 1. We are 
interested in the entanglement of a given initially pure 
state 1 4*) as a function of time. That is, the initial state 
suffers from decoherence described by this model and is 
given after a certain time t by 



p(t) =£l£2...fAr|*)(*| 



(2) 



We will use as a criterion the distillable entanglement 
of the system, which tells one whether it is possible to cre- 
ate true (irreducible) multiparticle entangled pure states. 
That is, we consider N distinct parties each holding a 
particle belonging to the TV-particle state p. The state 
p is called iV-party distillable entangled if there exists a 
local protocol (i.e. each of the parties acts independently 
on their systems) such that one can obtain from a suffi- 
ciently large number of copies of p some true TV-particle 
entangled pure state f^. Note that in this definition it 
is not necessary to specify the kind of multiparticle en- 
tangled pure state which is created, as all true TV-party 
entangled pure states can be obtained from each other if 
several copies of the state are available We remark 
that a necessary condition for TV-party distillability is 
that the partial transpositions with respect to any group 
of parties are non-positive :4j . 

We start out by investigating the entanglement prop- 
erties of GHZ states under this decoherence model. 
One readily finds that \GHZ) evolves to a state 
p{t) diagonal in the GHZ-basis {l^^^^j^ few-i) ~ 

1/V2(|fcifc2 . . . fcAT-iO) ± \hh ■ ■ ■ ^jv-il))}, fcj e {0, 1}, 
with coefficients that depend on fc = kj . For k ^ 

we have = A^T = Afe, where 



A, = [(1 + pni - pY'-'^ + {I+pY'-^I - p)'=]/2^^+\(3) 

while A+ = Ao + j3"/2, Ag = Aq - p^/2. Note that 
Ai > A2 > . . . > A[Ar/2]. The partial transposition with 



respect to a group which contains exactly k parties is 
positive, p{tY^k > 0, if and only if 3 



P 



N 



< 2Xk 



(4) 



For k = 1 one observes that the threshold value Pcrit = 
g-Ktcrit -vvhere the partial transposition with respect to 
one particle becomes positive increases with TV. This im- 
plies that for t > tcrit = T the state is no longer TV-party 
distillable entangled and thus the lifetime r of true TV- 
party entanglement decreases with the size of the system 
as expected (see Fig. 1). 
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FIG. 1: Fig. 1(a): Upper bound on lifetime kt of Af-party 
entanglement in systems with TV ^ oo particles for different 
M. Fig. 1(b): Same as Fig. 1(a) but with double-logarithmic 
axis. Note that the same figures are obtained for the lifetime 
of iV-particle entanglement in TV-particle systems (see text). 

We now turn our attention to a class of multiparti- 
cle entangled states, the so called graph states (fl]. This 
class includes a variety of entangled states, e.g. GHZ 
states, cluster states and codewords of error correction 
codes Graph states, l^^ji/js-.-ziw)! ^^'^ defined as joint 
eigenstates of a set of TV commuting correlation operators 
Kj associated with a graph G = {V, E) which is a set of 
TV vertices V connected in a specific way by edges E. The 
edges specify the neighborhood relation between vertices, 

and Kj are given by Kj = Ux'' ]\{k,j}£E '^^^^ ■ The graph 
states \^ ^l^^l^...^lfY) fulfill the set of eigenvalue equations 
/^,|«'miP2...p«> - (-ir'l*MiP2...p«> Vj, p, e {0,1} and 
form a basis in H = ((T^)®^. 

In the following, we consider a linear cluster state of 
TV qubits @ specified by the graph G with edges (A;, A; -|- 
l)Vfc and investigate the entanglement properties of this 
state under the decoherence model described by Eq. ^ . 
We will show that the lifetime of distillable TV-particle 
entanglement is for this state independent of the size of 
the system. In particular, we establish a lower bound t< 
on the lifetime iciit of distillable TV-party entanglement 
for such states which is independent of the number of 
particles TV. This is in sharp contrast to the behavior of 
GHZ states. 

In order to distill TV-party entangled states, it is suffi- 
cient that maximally entangled pairs between all neigh- 
boring parties can be created, which we will show in the 
following. We emphasize that we use the distillability of 
neighboring pairs only as a tool to prove TV-party distil- 
lability and it does not mean that the entanglement of 
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the cluster state were in some sense only "bipartite" . We 
make use of the following properties of graph states: (i) 
Measuring all but two neighboring particles, say fc, j of 
a graph state in the eigenbasis of cr^ results in the 
creation of another graph state with only a single edge 
{k,j} 8]. The resulting state of particles k,j is up to 
local az operations equivalent to a maximally entangled 
state of the form |$) ee 1/\/2(|0)^|0)^ + , where 

\k)x [\k)z] denote eigenstates of ax [(Jz] respectively, (ii) 
The action of a depolarizing channel £k on a graph state 
can equivalently be described by a map M k whose Kraus 
operators only contain products of Pauli matrices CTz and 
the identity, where here CTz may act on particles k and 
its neighbors, i.e. particles which are (in the correspond- 
ing graph) connected by edges to particle k. This follows 
from the fact that fii^^l^^r) = {-lY' ai^^ Kjl-^i p) , where 

ax Kj is an operator which contains only products of cr^ 
operators at neighboring particles of particle j, and the 
identity otherwise. Similarly, the action of dy"'^ on graph 
states is up to a phase factor equivalent to the action of 
an operator which contains only products of cr^ operators 
acting on particle j and all its neighbors. 

A sufficient condition when bipartite entanglement be- 
tween two neighboring particles, say particle k and fc -I- 1, 
can be created from p{t) can be found straightforwardly. 
To this aim, one performs measurements in the eigenbasis 
of az on all but particles k and k+1 (We remark that mea- 
surements on all neighboring particles of particles k,k+l 
would also be sufficient) . From properties (i) and (ii) fol- 
lows that these measurements commute with the action 
of the CPM on the cluster state. The resulting state can 
thus be described by M1M2 ■ ■ ■ MN\^)k,k+i{^\'S)\x){x\, 
where |x) is a state of the remaining (A^— 2) particles, and 
1$) is a maximally entangled state equivalent up to az 
operations to |$). It is important to note that A4j acts 
trivially on particles k and k+1 if j ^ {fc— 1, k, k+1, k+2} 
since any kind of errors in graph states only effect the 
corresponding particles and/or its neighbors (see (ii)). 
This implies that the resulting state pk,k+i of particles 
k and k + 1 after tracing out the remaining particles 
can be obtained by considering only the (reduced) action 
of Mk-i,Mk,Mk+i,Mk+2 on the state |$)($|. One 
has that pk,k±\ is distillable if its partial transpose is 
non-positive [3 and finds a threshold value p< = 0.717 
(k<< = 0.3327). We emphasize that this threshold value 
is independent of iV, as only errors acting on the direct 
neighborhood of the pair of particles in question influ- 
ences the threshold value. Since this method allows one 
to distill maximally entangled pairs between arbitrary 
pairs of neighboring particles, we have that for p > P< 
the linear cluster state remains A^~party distillable, in- 
dependent of the number of particles iV. 

We point out that above results are not restricted to 
linear cluster states but similarly hold for all graph states 
associated with some lattice geometry. To be specific, if 
we consider a family of graph states whose maximum de- 
gree does not increase with the number of particles N and 



two neighboring vertices k \j\ of the corresponding graph 
with disjoint sets oiuk + l [uj + 1] neighbors respectively, 
one can distill a maximally entangled pair between the 
two parties if pV4[(l -Hp"'=)(l + (1 -P^)/4 > 1/2. 
A similar expression can be obtained for more general 
graphs where k and / may have common neighbors. It 
turns out that ii m = dk + dj — 2, where di is the degree 
of vertex i, then for nt < — ln(i)/([Y] + 2) the state 
is certainly distillable. For cluster states correspond- 
ing to a regular 2D [3D] lattice one finds p< — 0.8281 
(Kt< = 0.1886) [p< = 0.8765 (Kt< = 0.1318)] respec- 
tively. This dependence on the degree is consistent with 
the results obtained for GHZ states, which correspond to 
a graph with edges (l,/c)Vfc. The degree of this graph 
increases with A^ and thus the corresponding lifetime de- 
creases. However, we emphasize that the dependence on 
the degree of the graph only applies to the lower bound t< 
established with this specific purification method and the 
general dependence of the lifetime of an arbitrary graph 
state on the (degree of the) graph is presently unknown. 

We remark that the observed behavior, i.e. that the 
lifetime of multiparticle entanglement for cluster (and 
similar) states is independent of the size of the system, 
also holds for more general decoherence models. This 
follows from the fact that — similar to (ii) — the action 
of any CPM acting on graph states describing an arbi- 
trary decoherence process can be estimated by a CPM 
whose Kraus operators only contain products of az oper- 
ators and the identity and thus the measurement (i) still 
commutes with the CPM. To this aim, we apply after 
the application of the CPM a local depolarization pro- 
cedure which maps arbitrary density operators to oper- 
ators diagonal in the graph state basis without changing 
the diagonal elements [Tj. When restricted on graph 
states, the resulting action of the initial CPM (given by 
J2k I '^k.iOkpOi, where Ok,Oi are products of Pauli op- 
erators)) can then be described by a CPM specified by 
J2k '^k.kOkpOj., where all operators Ok can be expressed 
in terms of products of az operators. Only operators 
Ok which act non-trivially on particles k,k + 1 or their 
neighbors in the graph affect the resulting maximally en- 
tangled pair after the measurement (i), leading again to 
a threshold value which is independent on the size of the 
system for all those decoherence models where the num- 
ber of such operators Ok is independent of A^. This is for 
instance the case if each Ok acts non-trivially on a finite, 
localized number of subsystems. 

In the following we investigate the effective size of the 
system that remains entangled during the decoherence 
process. The fact that A^-party entanglement vanishes 
after a certain time does not imply that all entangle- 
ment has disappeared. We consider partitions of the A^- 
particle state into M groups, where now parties within 
one group are allowed to perform joint operations. For 
states associated to some lattice geometry, specific par- 
titions correspond to a re-scaling of the size of the sub- 
systems. A A^-particle state is called Af -party distillable 
entangled if there exists some partitioning into M par- 
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ties (M-partitioning) such that the state is distillable to 
some true M-party entangled pure state Q . If the state 
is not A^-party distillable entangled, it can well be that 
it is M-party distillable entangled for some M < N. We 
will be interested in the lifetime of M-party distillable 
entanglement. 

For GHZ states, one can indeed determine analytically 
the lifetime of M-party entanglement for any M. Recall 
the condition for positivity of partial transposition with 
respect to a group of k parties given by Eq. Since 
> Aj for j < i, we have that the group containing the 
fewest number of parties determines the threshold value 
for which the state is no longer M-party distillable, as the 
corresponding partial transposition with respect to this 
group is the first one to become positive. Thus M-party 
entanglement corresponding to a specific M-partitioning 
has longest lifetime if all groups have (approximately) the 
same size. For a minimal group size of m particles, a TV- 
particle GHZ state can contain at most M = [N/ni] such 
groups of size m. This allows one to obtain the maximum 
lifetime of M-party entanglement which is determined by 
Eq. with fc = TO. 

One obtains an upper bound on the lifetime of M- 
party entanglement if one approximates Xm by some 
A™ < Xm and chooses A™ = (1 - +p)^""V2^+i. 
In this case, the partial transposition with respect to a 
group of TO parties is certainly positive if 2A,„ > , 
which can be rewritten asm < N log[2p/ (1 +p)]/ log[(l — 
p)/{l+p)]. Using that M = [N/m] , we have that a iV- 
particle state is certainly no longer M-party entangled 
if 

M > [log(l -p) - log(l +p)]/[log(2p) - log(l (5) 

On the one hand, Eq. lO — illustrated in Fig. 1 — pro- 
vides an upper bound on the lifetime ktm of M-party 
entanglement in the system. On the other hand, for a 
fixed time t Eq. Q allows one to determine the maxi- 
mum M of distillable multiparty entanglement remain- 
ing in the system. One observes (see Fig. 1) that M 
rapidly decreases with t. For nt ^ 1, one finds that 
M « -21og(Ki)/(Ki), while for Kt > 0.8049, we have 
that also 2-party entanglement disappears and the state 
becomes fully separable as all partial transposes are pos- 
itive (which is a sufficient condition for separability for 
states diagonal in the GHZ state basis 4] ) . We also point 
out that this bound on the maximal size of M-party en- 
tanglement is independent on the number of particles N. 
That is, even if one considers groups of size to — > oo 
(for a total number of particles N — > cxd), the maximum 
number M of such groups which remain entangled after 
a certain time t is finite. While for the finest partition 
of the system into N parties the lifetime ktn essentially 
scales like 1/N, for the coarsest partition of the N par- 
ties into two groups we have no scaling behavior with N, 
i.e. KT2 ~ const. It follows that in the limit N —^ oo any 
partitioning in groups of finite size m leads to a vanishing 
lifetime of the corresponding M = N/ m party entangle- 
ment. Only if one considers a fixed number of groups 



M whose size to = N/M grows to infinity as iV ^ oo, 
i.e. the groups itself are of macroscopic size, one obtains 
a finite lifetime of the corresponding M-party entangle- 
ment. 

We remark that Eq. ||SJ| also allows one to deter- 
mine the lifetime of genuine iV-party entanglement, ob- 
tained from M = N/m by setting to = 1. In a sim- 
ilar way, one can obtain a lower bound on the life- 
time of M-party entanglement by choosing A^ = 2(1 — 
p)™(l -f p)^-™/2^+i > X,n- One finds that for M < 
log[2{l - p) / {1 + p)]/ \og[2p/ {1 + p)] all partial transpo- 
sition with respect to this M-partitioning are certainly 
non-positive, which for these kind of states already en- 
sures that the state is M-party distillable . The condi- 
tion for 2-party entanglement has recently been derived 
by Simon and Kempe in Ref . , who observed that the 
threshold value for p such that the partial transposition 
with respect to a partition N/2-N/2 parties is positive, 
decreases with the size of the system. Based on this ob- 
servation, they conclude that GHZ states of more parti- 
cles are more stable against local decoherence. However, 
as pointed out in this paper, the effective number of sub- 
systems that remain entangled decreases with time, such 
that the entanglement becomes bipartite when approach- 
ing the threshold value found by Simon and Kempe. The 
lifetime of genuine iV-party entanglement thus decreases 
in fact with the size of the system. 

On the other hand, for cluster states (and similar graph 
states) one can show that there is no scaling with re- 
spect to either the size of the partitions or N. In this 
sense both the kind of entanglement in the system and 
its lifetime are universal. To this aim, it is sufficient to 
show that there exists not only a lower bound but also 
an upper bound on the lifetime of entanglement which is 
independent of N and valid for any partition. To be spe- 
cific, we show that for any cluster state there exist times 
t>, i< independent of N such that for t > the state is 
separable with respect to the finest partition (and hence 
not distillable with respect to any partition), while for 
t < i< the state is distillable with respect to the finest 
partition (and hence distillable with respect to any par- 
tition). The existence of such a time t< was already 
shown earlier in this paper, so we concentrate on i>. 
We use that (i) the graph state l^'g) corresponding to a 
graph G can be written as (si l^-) = U{k,i)eE Uki\+)'^'^ , 
where Uki = d'^'+^f )/2®(l'"--<")/2 and |+) - 
1/V2(|0> + |1>); (u) The map £k{p) (Eq. ^ with noise 
parameter p = e"*^') can equivalently be described by 

a concatenation of three maps £]^\px)£^\py)£'l^\pz) 

where £i'\pj)p = PjP + (1 ~ p,)/2[p + af^paf^] 
with pj ~ J = x,y,z. We consider the map 

£ii\p,)p = £i'\p,)£['\p,)UupUl. Using the re- 

suits of Ref. HU], one can verify that £j.i {pz) is 
separable and hence not able to create entanglement 
if Pz < (\/2 — 1). Consider now a graph state 
where each vertex in the corresponding graph has to 
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neighbors. Since Uki commutes with £j (pz) and 
^k\p^)P = ^k\VPd^k'\\/P^)P^ "^e can rewrite p{t) = 

nf=ifr(p.)ff (p.)n(.,).^flf(pi^/™)i+)(+i«^. We 

have that all maps in this expression are separable if 
Pz^™ < (\/2 — 1) which implies that p{t) is fully separable 
for Kt > — 2m ln(-\/2 — !)• This provides the announced 
upper bound t> on the lifetime of distillable entangle- 
ment valid for arbitrary partitions. 

In this paper, we have studied the behavior of multi- 
particle entangled states under decoherence. For GHZ 
states, we found that the lifetime of true M-party en- 



tanglement decreases with the size of the system and the 
effective number of entangled subsystems decreases with 
time. For cluster and similar graph states, however, we 
have shown that the lifetime of TV-party entanglement 
is independent of the size of the system. These results 
suggests that true multiparticle entanglement in macro- 
scopic objects can be more stable (and might be more 
common) than previously thought of. 

We thank J.I. Cirac for valuable discussions. This work 
was supported by European Union (HPMF-CT-2001- 
01209 (W.D.), IST-2001-38877,-39227) and the DFG. 



[1] E. Schrodinger, Die Naturwissenschaften 23, 807 (1935). 
[2] For a recent review see e.g. W. Zurek, quant-ph/0306072 
[3] Fortschritte der Physik, Vol. 48, Number 9-11 (2000). ; 
[4] W. Diir and J. I. Cirac, Phys. Rev. A 61, 042314 (2000); 

W. Diir and J. I. Cirac, Phys. Rev. A 62, 022302 (2000). 
[5] A. V. Thapliyal and J. A. Smolin, quant-ph/0212098 
[6] H.-J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 

910 (2001). R. Raussendorf, D. Browne and H.-J. Briegel, 

quant-ph/0301052 
[7] C. Simon and J. Kempe, Phys. Rev. A 65, 052327 (2002). 



[8] M. Hein, J. Eisert and H.-J. Briegel, quant-ph/0307130| 
[9] M. Horodecki, P. Horodecki and R. Horodecki, Phys. 

Rev. Lett. 78, 574 (1997). 
[10] D. Schlingemann and R. F. Werner, Phys. Rev. A 65, 

012308 (2002). 

[11] W. Diir, H. Aschauer and H.-J. Briegel, 

quant-ph/0303087 
[12] J. I. Cirac et al. Phys. Rev. Lett. 86, 544 (2001); 



